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(54) CONTROL DEVICE OF INTERNAL COMBUSTION ENGINE 



(57) A control system of an internal combustion en- 
gine provided with a throttle valve passage air flow cal- 
culation equation by which the throttle valve passage air 
flow mt is expressed as a function of the downstream 

side intake pipe pressure at the downstream side of the 
throttle valve and a cylinder intake air flow calculation 
equation by which the cylinder intake air flow mc is ex- 
pressed as a function of the downstream side intake pipe 



pressure, where said downstream side intake pipe pres- 
sure Pm and cylinder intake air flow mc when the throttle 
valve passage air flow mt found from said throttle valve 
passage air flow calculation equation and the cylinder 
intake air flow mc found from said cylinder intake airflow 
calculation equation match are calculated as the down- 
stream side intake pipe pressure Pmta and cylinder in- 
take air flow mcta at the time of steady operation under 
the operating conditions at that time is provided. 
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Description 
TECHNICAL FIELD 

5 [0001] The present invention relates to a control system of an internal connbustlon engine. 
BACKGROUND ART 

[0002] In recent years, studies have been conducted on constructing a nnodel of the intake system of an internal 

10 combustion engine leased on fluid dynamics etc. and using that model to calculate parameters to control the internal 
combustion engine. That is, for example, a throttle model, intake pipe model, intake valve model, etc. have been con- 
structed for the intake system of internal combustion engines, these models have been used to calculate the cylinder 
filling air amount etc. from the throttle valve opening degree, atmospherk: pressure, atmospheric temperature, etc. and 
these have been used to control internal combustion engines. 

15 [0003] However, when controlling an intemal combustion engine, in particular when using the above model to control 
an intemal combustion engine, to calculate the parameters related to the control, sometimes the throttle valve downstream 
side intal<e pipe pressure Pmta or cylinder Intake airflow mcta at the time of steady operation (or the cylinder air filling 
rate KIta at the time of steady operation obtained calculated from that (that is, the ratio of mass of cylinder filling air with 
respect to the mass of the air of the total stroke capacity of the cylinder)) becomes necessary. For example, Japanese 

20 Patent Publication (A) No. 2001-41 095 discloses a method of calculating a throttle valve passage air flow based on the 
throttle valve downstream side intake pipe pressure, atmospheric pressure, etc. at that time and the Pmta. 
[0004] Further, the above throttle valve downstream side intake pipe pressure Pmta or cylinder intake air flow mcta 
at the time of steady operation has conventionally been found using a map. That is, for example, in the Japanese Patent 
Publication (A) No. 2001-41095, the Pmta is found from a map using the throttle valve opening degree and/or engine 

25 speed etc. as arguments. 

[0005] However, when actually preparing such a map, a tremendous amount of time becomes necessary. That is, to 
prepare a map, it is necessary to actually measure the Pmta or mcta while successively changing the arguments. This 
work becomes tremendous. Further, there is the concern that an increase in the necessary maps or arguments will 
increase the map searching operation and increase the control load. 

30 

DISCLOSURE OF THE INVENTION 

[0006] The present invention was made in consideration of this problem and has as its object the provision of a control 
system of an internal combustion engine designed to find at least one of the throttle valve downstream side intake pipe 
35 pressure Pmta and cylinder intake air flow mcta at the time of steady operation by a simpler method. 

[0007] The present invention provides as a means for solving the above problem a control system of an internal 
combustion engine as described in the claims of the claim section. 

[0008] In a first aspect of the present invention, there is provided a control system of an internal combustion engine 
provided with a throttle valve passage airflow calculation equation by which a throttle valve passage air flow is expressed 

<o as a function of a downstream side intake pipe pressure at the downstream side of a throttle valve and a cylinder intake 
air flow calculation equation by which a cylinder intake airflow is expressed as a function of the downstream side Intake 
pipe pressure, the downstream side intake pipe pressure when the throttle valve passage air flow found from the throttle 
valve passage airflow calculation equation and the cylinderintake air flow found from the cylinder intake air flow calculation 
equation match being calculated as the downstream side intake pipe pressure at the time of steady operation under the 

45 operating conditions at that time. 

[0009] The downstream side intake pipe pressure at the time of steady operation had in the past been found using a 
map, but there were the problems that the manhours involved in mapmaking work was tremendous and the control load 
at the time of map searching was also large. 

[0010] As opposed to this, in this aspect, the fact that at the time of steady operation, the throttle valve passage air 
50 flow and the cylinder intake air flow match is utilized to find by calculation the downstream side intake pipe pressure at 
the time of the steady operation. For this reason, according to the present embodiment, it is possible to more simply find 
the downstream side intake pipe pressure at the time of steady operation. 

[0011] In a second aspect of the present invention, there is provided a control system of an internal combustion engine 
provided with a throttle valve passage airflow calculation equation by which a throttle valve passage airflow is expressed 
55 as a function of a downstream side intake pipe pressure at the downstream side of a throttle valve and a cylinder intake 
air flow calculation equation by which a cylinder intake air flow is expressed as a function of the downstream side intake 
pipe pressure, the cylinder Intake air flow when the throttle valve passage air flow found from the throttle vatve passage 
air flow calculation equation and the cylinder intake air flow found from the cylinder intake air flow calculation equation 
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match being calculated as the cylinder intake air flow at the time of steady operation under the operating conditions at 
that time. 

[0012] The cylinder intake air flow at the time of steady operation had in the past been found using a map. There was 
the same problem when finding the above downstream side intake pipe pressure at the time of steady operation by a map. 
[0013] As opposed to this, in this aspect, the fact that at the time of steady operation, the throttle valve passage air 
flow and the cylinder intake air flow match is utilized to find by calculation the cylinder intake air flow at the time of the 
steady operation. For this reason, according to the present embodiment, it is possible to more simply find the cylinder 
intake air flow at the time of steady operation. 

[0014] In a third aspect of the present invention, the cylinder intake air flow when a throttle valve passage air flow 
found from the throttle valve passage air flow cateulation equation and a cylinder intake airflow found from the cylinder 
intake air flow calculation equation match is calculated as the cylinder intake air flow at the time of steady operation 
under the operating conditions at that time. 

[0015] According to the present embodiment, it is possible to more simply find the downstream side intake pipe 
pressure and/or cylinder intake air flow at the time of steady operation. 

[001 6] In a fourth aspect of the present Invention, the throttle valve passage air flow calculation equation is expressed 
as the following equation (1 ) where mt is a throttle valve passage air flow, jji is a flow coefficient at the throttle valve, At 
is a cross-sectional area of the opening of the throttle valve, Pa is an atmospheric pressure, Ta is an atmospheric 
temperature, R is a gas constant, Pm is a downstream side intake pipe pressure, and 0(Pm/Pa) is a coefficient determined 
in accordance with the value of Pm/Pa, and the cylinder Intake air flow calculation equation is expressed as the following 
equation (2) where mc is the cylinder intake air flow and a and b are compliance parameters detemnined based on at 
least the engine speed: 



[0017] According to the present embodiment, relatively simple calculation may be used to accurately find the down- 
stream side intake pipe pressure and/or cylinder intake air flow at the time of steady operation. 
[0018] In a fifth aspect of the invention, the internal combustion engine has an exhaust gas recirculation passage for 
making at least part of the exhaust gas discharged into the exhaust passage flow into the intake passage and an EGR 
control valve for adjusting the flow of the exhaust gas passing through the exhaust gas recirculation passage, the throttle 
valve passage air flow calculation equation is expressed as the following equation (3) wherein mt is a throttle valve 
passage air flow, jx is a flow coefficient at the throttle valve, At is a cross-sectional area of the opening of the throttle 
valve, Pa is an atmospheric pressure, Ta is an atmospheric temperature, R is a gas constant, Pm is a downstream side 
intake pipe pressure, and <l>(Pm/Pa) is a coefficient determined in accordance with the value of Pm/Pa, and the cylinder 
intake air flow calculation equation is expressed as the following equation (4) where mc is a cylinder intake air flow, and 
e and g are compliance parameters detennined based on at least an engine speed and an opening degree of the EGR 
control valve. 



[0019] According to the present embodiment, even when perfonning exhaust gas recirculation, a relatively simple 
calculation may be used to accurately find the downstream side intake pipe pressure and/or cylinder intake air flow at 
the time of steady operation. 





mc = e • Pm+g 



(4) 
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[0020] In a sixth aspect of the invention, the Internal combustion engine further has a variable valve timing mechanism 
for changing an operating timing of a valve provided in each cylinder and, based on the compliance parameters e and 
9 when the operating timing is a first valve timing and the EGR control vaive is at a first opening degree, the compliance 
parameters e and g when the operating timing is the first valve timing and the EGR control valve is at a second opening 
5 degree, and the compliance parameters e and g when the operating timing is a second valve timing and the EGR control 
valve is at a first opening degree, the compliance parameters e and £ when the operating timing is the second valve 
timing and the EGR control valve is at the second opening degree are estimated. 

[0021] According to the present embodiment, when exhaust gas recirculation is perfomried and there is a variable 
valve timing mechanism, it is possible to reduce the manhours of the mapmaking work for the compliance parameters 
10 e and fl. Further, if reducing the number of stored maps, it is also possible to reduce the control load at the time of map 
searching. 

[0022] In a seventh aspect of the present invention, when the compliance parameters e and 3 when the operating 
timing is the second valve timing and the EGR control valve is at a first opening degree respectively take two different 
values when the throttle valve downstream side intake pipe pressure is larger than and smaller than a first pressure and 

15 the compliance parameters e and 3 when the operating timing is the second valve timing and the EGR control valve is 
at the second opening degree are estimated to take three or more different values in accordance with the throttle valve 
downstream side intake pipe pressure, based on the compliance parameters e and g when the operating timing is a first 
valve timing and the EGR control valve is at a first opening degree, the compliance parameters e and g[ when the 
operating timing is the first valve timing and the EGR control valve is at a second opening degree, and the compliance 

20 parameters e and £ when the operating timing is a second valve timing and the EGR control valve is at a first opening 
degree, approximated compliance parameters ep and gp designed to take two values differing when the throttle valve 
downstream side intake pipe pressure is larger and smaller than a first pressure are calculated and these are made the 
compliance parameters e and c[ when the operating timing Is the second valve timing and the EGR control valve is at 
the second opening degree. 

25 [0023] According to the present embodiment, the processing when finding the downstream side intake pipe pressure 

and/or cylinder intake air flow at the time of steady operation is simplified and the control load can be reduced. 
[0024] In an eighth aspect of the invention, the case where the EGR control valve is at the first opening degree is the 
case where the EGR control valve is closed. 

[0025] By using the case where the EG R control valve is closed as a standard, it is possible to more accurately estimate 
30 the compliance parameters e and g in the case where the operating timing is the second (that is, any) valve timing and 
the EGR control valve is at the second (that is, any) opening degree. Further, as a result, it is possible to more accurately 
find the downstream side intake pipe pressure and/or cylinder intake air flow at the time of steady operation. 
[0026] In a ninth aspect of the invention, at a portion where the throttle valve passage air flow mt and cylinder intake 
air flow mc invert in magnitude, the throttle valve passage air flow calculation equation used is an approximation equation 
35 expressed by a linear equation of the downstream side intake pipe pressure Pm. 

[0027] In a 10th aspect of the invention, the approximation equation is made a linear equation expressing a line 
connecting two points on a curve expressed by the throttle valve passage airflow calculation equation and before and 
after the point where the throttle valve passage air flow mt and cylinder intake air flow mc invert in magnitude. 
[0028] According to the ninth and 1 0th aspects, the calculation when finding the downstream side intake pipe pressure 
^ and/or cylinder intake air flow at the time of steady operation is simplified and the control load can be reduced. 

[0029] In an 11 th aspect of the invention, instead of the atmospheric pressure Pa, a throttle valve upstream side intake 
pipe pressure Pac found considering at least a pressure loss of an air cleaner is used. 

[0030] According to the present embodiment, it is possible to more accurately find the downstream side intake pipe 
pressure and/or cylinder intake air flow at the time of steady operation. 

^5 [0031] In a 12th aspect of the invention, a throttle valve upstream side intake pipe pressure Pac found considering at 
least a pressure loss of an air cleaner is found based on the previously found throttle valve passage air flow, and the 
approximation equation is made a linear equation expressing a line connecting two points shown by coordinates obtained 
by multiplying Pac/Pa with values of the downstream side intake pipe pressure and the throttle valve passage airflow 
showing coordinates of two points on a curve expressed by the throttle valve passage airflow calculation equation and 

so before and after a point where a throttle valve passage air flow mt and cylinder intake air flow mc invert in magnitude. 
[0032] According to the present embodiment, the calculation when finding the downstream side intake pipe pressure 
and/or cylinder intake air flow at the time of steady operation is simplified and the control load can be reduced. Further, 
by considering the pressure loss of the air cleaner etc., it is possible to more accurately find the downstream side intake 
pipe pressure and/or cylinder intake air flow at the time of steady operation. 

55 [0033] Below, the present invention will be understood more fully from the attached drawings and the description of 
the preferred embodiments of the present invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0034] 

5 FIG. 1 is a schematic view of an example of the case of applying a control system of an internal combustion engine 

of the present Invention to a cylinder injection type spark ignition internal combustion engine. 
FIG, 2 is a view of an intal<e air amount model. 

FIG. 3 is a view of the relationship between the throttle valve opening degree and flow coefficient. 
FIG. 4 is a view of a function «>(Pm/Pa). 
10 FIG. 5 is a view of the basic concept of a throttle model. 

FIG. 6 is a view of the basic concept of an intake pipe model. 
FIG. 7 is a view of the basic concept of an intake valve model. 

FIG. 8 is a view of definitions of the cylinder filling air amount and cylinder intake airflow. 
FIG. 9 is a view of the relationship of a downstream side intake pipe pressure Pm and throttle valve passage air 
15 flow mt and cylinder intake air flow mc and shows that the downstream side intake pipe pressure Pm when the 

throttle valve passage air flow mt and cylinder intake air flow mc become equal is the downstream side intake pipe 
pressure Pmta at the time of steady operation and that the cylinder intake airflow mc at that time is the cylinder 
intake air flow mcta at the time of steady operation. 

FIG. 1 0 is an enlarged view of the vicinity of an intersecting point EP for a view similar to FIG. 9 and is a view to 
20 explain approximating a curve expressing the throttle valve passage air flow mt by a line and approximating two 

lines expressing the cylinder intake air flow mc by a single line. 

FIG. 11 is a schematic view of an example of the case of applying a control system of an internal combustion engine 
of the present invention to a cylinder injection type spark ignition internal combustion engine different from FIG. 1 . 
FIG. 12 is aviewforexplainingthe method of utilizing compliance parameters e and fl under predetemnined conditions 
25 to estimate the compliance parameters e and q under any conditions. 

FIG. 13 is also a view for explaining the method of utilizing compliance parameters e and under predetermined 
conditions to estimate the compliance parameters e and g under any conditions. 

FIG. 14 is a view for explaining the method of approximating a cylinder intake air flow mcl 1 expressed by three 
lines by an approximated cylinder intake air flow me'11 expressed by two lines and shows the case where the 
30 predetemiined pressure Pml is larger than a predetennined pressure Pm2. 

FIG. 15 is a view similar to FIG. 14 and shows the case where the predetennined pressure Pm1 Is smaller than a 
predetennined pressure Pm2. 

BEST MODE FOR WORKING THE INVENTION 

35 

[0035] Below, embodiments of the present Invention will be explained in detail with reference to the drawings. Note 
that in the figures, the same or similar elements are assigned common reference numerals. 

[0036] FIG. 1 is a schematic view showing an example of the case of applying a control system of an internal combustion 
engine of the present invention to a cylinder injection type spark ignition internal combustion engine. Note that the present 
^ invention may also be applied to another spark ignition type internal combustion engine or compression ignition type 
internal combustion engine. 

[0037] As shown in FIG. 1, an engine body 1 is provided with a cylinder block 2, pistons 3 reciprocating inside the 
cylinder block 2, and a cylinder head 4 fixed on the cylinder block 2. The pistons 3 and cylinder head 4 fomi combustion 
chambers 6 between them. The cylinder head 4 is provided with, for each cylinder, intake valves 6, intake ports 7, 
45 exhaust valves 8, and exhaust ports 9. Further, as shown in FIG. 1, sparkplug 10 is provided at the center of the inside 
wall surface of the cylinder head 4, while fuel injector 1 1 is provided at the periphery of the inside wall surface of the 
cylinder head 4. Further, the top surface of each piston 3 is fonned with a cavity 1 2 extending from below the fuel injector 
1 1 to below the spark plug 1 0. 

[0038] The intake ports 7 of each cylinder are connected through a downstream side intake pipe 13 to a surge tank 
50 1 4, while the surge tank 1 4 is connected through an upstream side intake pipe 1 5 to an air cleaner 1 6. Inside the intake 
pipe 1 5, a throttle valve 1 8 driven by a step motor 1 7 is provided. On the other hand, the exhaust ports 9 of each cylinder 
are connected to the exhaust pipe 19, and this exhaust pipe 19 is connected to an exhaust purifier 20. 
[0039] The electronic control unit (ECU) 31 is comprised of a digital computer provided with a RAM (random access 
memory) 33, ROM (read only memory) 34, CPU (microprocessor) 35, input port 36, and output port 37, connected with 
55 each other by a bidirectional bus 32. The intake pipe 1 3 is provided with an intake pipe pressure sensor 40 for detecting 
the pressure in the intake pipe. The intake pipe pressure sensor 40 generates an output voltage proportional to the 
intake pipe pressure. This output voltage is input through the corresponding AD converter 38 to the input port 38. 
[0040] Further, a throttle valve opening degree sensor 43 for detecting the opening degree of the throttle valve 1 8. an 
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atmospheric pressure sensor 44 for detecting the pressure of the atmosphere around the internal combustion engine 
or the pressure of the air taken into the intake pipe 1 5 (intake pressure), and an atmospheric temperature sensor 45 for 
detecting the temperature of the atmosphere around the intemal combustion engine or the temperature of the air taken 
into the intake pipe 1 5 (intake temperature) are provided. The output voltages of these sensors are input to an input port 
36 through corresponding AD converters 38. Further, the accelerator pedal 46 is connected to a load sensor47 generating 
an output voltage proportional to the amount of depression of an accelerator pedal 46, while the output voltage of the 
load sensor 47 is input to an input port 36 through the corresponding AD converter 38. A crank angle sensor 48 generates 
an output pulse each time for example the crankshaft rotates by 30 degrees. This output pulse is input to the input port 
36. The CPU 35 uses the output pulse of this crank angle sensor 48 to calculate the engine speed. On the other hand, 
the output port 37 is connected through a corresponding drive circuit 39 to the spark plugs 10, fuel injectors 11, step 
motor 1 7, etc. 

[0041] Note that, in recent years, control systems of internal combustion engines controlling intemal combustion 
engines based on parameters calculated using the models of the intake systems of internal combustion engines con- 
structed based on fluid dynamics, etc. have been studied. That is, for example, a throttle model, intake pipe model, 
intake valve model, etc. have been constructed for the intake systems of internal combustion engines, these models 
have been used to calculate the cylinder filling air amount etc. from the throttle valve opening degree, atmospheric 
pressure, atmospheric temperature, etc., and the internal combustion engine has been controlled based on this. 
[0042] Further, in the present embodiment as well, in the configuration shown in FIG. 1, a model is used to control 
the intemal combustion engine. That is, in the present embodiment, usually, the Intake air amount model M20 explained 
below Is used for control. FIG. 2 is a view of the intake air amount model M20. 

[0043] The intake air amount model M20, as shown in FIG. 2, is provided with a throttle model M21 , intake pipe model 
M22, and intake valve model M23. The throttle model M21 uses as input the opening degree et of the throttle valve 
detected by the throttle valve opening degree sensor (hereinafter referred to as the "throttle valve opening degree"), the 
atmospheric pressure Pa around the internal combustion engine detected by an atmospheric pressure sensor, the 
atmospheric temperature Ta around the internal combustion engine detected by an atmospheric temperature sensor, 
and the pressure Pm in the intake pipe at the downstream side from the throttle valve calculated in the later explained 
intake pipe model M22 (hereinafter referred to as the "downstream side intake pipe pressure"). By entering the values 
of these input parameters into the model equations of the later explained throttle model M21 , the flow of the air passing 
through the throttle valve per unit time (hereinafter referred to as the "throttle valve passage air flow mt") is calculated. 
The throttle valve passage air flow mt calculated in the throttle model M21 is input to the intake pipe model M22. 
[0044] The intake pipe model M22 uses as input the throttle valve passage air flow mt calculated in the throttle model 
M21 and the flow of air flowing into the combustion chamber per unit time explained In detail below (hereinafter referred 
to as the "cylinder intake air flow mc". Note that the definition of the cylinder intake air flow mc is described in detail in 
the intake valve model M23). By entering the values of these input parameters into the model equations of the later 
explained intake pipe model M22, the downstream side Intake pipe pressure Pm and the temperature Tm in the intake 
pipe at the downstream side of the throttle valve (hereinafter referred to as the "downstream side intake pipe temperature") 
are calculated. The downstream side intake pipe pressure Pm calculated at the intake pipe model M22 is input to the 
intake valve model M23 and throttle model M21 . 

[0045] The intake valve model M23 uses as input the downstream side intake pipe pressure Pm calculated at the 
intake pipe model M22. By entering this value into the model equations of the later explained intake valve model M23, 
the cylinder intake air flow mc is calculated. The calculated cylinder intake air flow mc is converted to the cylinder filling 
air amount Mc. Based on this cylinder filling air amount Mc, the amount of fuel injection from the fuel injector is detemnined. 
Further, the cylinder intake air flow mc calculated at the intake valve model M23 is input to the intake pipe model M22. 
[0046] As will be understood from FIG. 2, in the intake air amount model M20, the value of the parameters calculated 
in a certain model are utilized as input values to another model, so in the intake air amount model M20 as a whole, the 
actually Input values are the throttle valve opening degree et, atmospheric pressure Pa. and atmospheric temperature 
Ta, that is, only three parameters. The cylinder filling air amount Mc is calculated from these three parameters. 
[0047] Next, each of the models M21 to M23 of the intake air amount model M20 will be explained. 
[0048] In the throttle model M21 , the throttle valve passage airflow mt(g/s) is calculated based on the following equation 
(5) from the atmospheric pressure Pa(kPa), atmospheric temperature Ta(K), downstream side intake pipe pressure Pm 
(kPa), and throttle valve opening degree et. Here, the |x in equation (5) is the flow coefficient in the throttle valve, is a 
function of the throttle valve opening degree Bt, and is determined from the map shown in FIG. 3. Further, At(m2) shows 
the cross-sectional area of the opening of the throttle valve (hereinafter referred to as "throttle opening area") and is a 
function of the throttle valve opening degree et Note that it is also possible to find fi-At combining these flow coefficient 
\i and throttle opening area At from the throttle valve opening degree 6t by a single map. Further, R is the gas constant. 
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mt=/i • At 



R • Ta \PaJ 



(5) 



70 



[0049] <D(Pm/Pa) is a function of the following equation (6). The k in this equation (6) is the specific heat ratio (ic=Cp 
(constant pressure specific heat)/Cv (constant volume specific heat), made constant value). This function <I>(Pnri/Pa) can 
be expressed by the graph as shown in FIG. 4, so this graph may be stored as a map in the ROM of the ECU and the 
value of <D(Pm/Pa) may be found from the map instead of actual calculation using equation (6). 
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[0050] Equation (5) and equation (6) of the throttle model M21 are obtained by making the pressure of the gas upstream 
of the throttle valve 18 the atmospheric pressure Pa, making the temperature of the gas upstream of the throttle valve 
1 8 the atmospheric temperature Ta, making the pressure of the gas passing through the throttle valve 1 8 the downstream 
side intake pipe pressure Pm, applying the Law of the Conservation of Mass, the Law of the Conservation of Energy, 
and the Law of Conservation of Motion to the model of the throttle valve 18 shown in FIG. 5, and utilizing the gas state 
equation, the definition of specific heat ratio, and Mayer's fomiula. 

[0051 ] In the intake pipe model M22, the downstream side intake pipe pressure Pm(kPa) and downstream side intake 
pipe temperature Tm(K) are calculated from the throttle valve passage air flow mt(g/s), cylinder intake air flow mc(g/s), 
and atmospheric temperature Ta(K) based on the following equation (7)and equation (8). Note that Vm(m3) in equation 
(7) and equation (8) is a constant equal to the volume of the portion 13' of the intake pipe etc. from the throttle valve to 
the intake valve (hereinafter referred to as the "Intake pipe part"). 



35 



— — = — • (rat-mc) 
dt^lmj Vm 



(7) 



40 



dPm R / ^ - T \ 

-K • — • (mt • Ta-mc • Tm) 

dt Vm 



(8) 



45 



SO 



[0052] Here, the intake pipe model M22 will be explained with reference to FIG. 6. If the total amount of gas of the 
intake pipe part 13' is M, the change over time of the total amount M of gas is equal to the difference between the flow 
of gas flowing into the intake pipe part 1 3', that is, the throttle valve passage air flow nnt, and the flow of gas flowing out 
from the intake pipe part 1 3', that is, the cylinder intake airflow mc. From the Law of Conservation of Mass, the following 
equation (9) is obtained. From this equation (9) and the gas state equation (Pm-Vm=M-R Tm), equation (7) is obtained. 



55 



dM , 
— = mt-mc 
dt 



(9) 



[0053] Further, the amount of change over time of the energy M-Cv-Tm of the gas of the intake pipe part 1 3' is equal 
to the difference between the energy of the gas flowing into the intake pipe part 13' and the energy of the gas flowing 
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out of the intake pipe part 13'. For this reason, if the temperature of the gas flowing into the intake pipe part 13* is made 
the atmospheric temperature Ta and the temperature of the gas flowing out from the intake pipe part 13' is made the 
downstream side intake pipe temperature Tm, aocording to the Law of the Conservation of Energy, the following equation 
(10) is obtained. From this equation (10) and the gas state equation, equation (8) is obtained: 



= Cp • mt • Ta-Cp • mc • Tm (1^) 

dt 



[0054] In the intake valve model M23, the cylinder intake air flow mc is calculated from the downstream side intake 
pipe pressure Pm based on the following equation (11). Note that a and b in equation (1 1) are compliance parameters 
detemiined based on at least the engine speed NE. A map is prepared in advance and the map is searched through in 
15 accordance with need to find these. 

mc=a • Pm-b (11) 

^ [0055] The abovennentioned intake valve model M23 will be explained with reference to FIG. 7. In general, the amount 
of air filled in the combustion chamber 5 when the intake valve 6 is closed, that is, the cylinder filling air amount Mc, is 
detennined when the intake valve 6 is closed (at the time the intake valve is closed) and is proportional to the pressure 
in the combustion chamber 5 at the time the intake valve is closed. Further, the pressure in the combustion chamber 5 
at the time the intake valve is closed can be deemed equal to the pressure of the gas upstream of the intake valve, that 

25 is, the downstream side intake pipe pressure Pm. Therefore, the cylinder filling air amount Mc can be approximated as 
being proportional to the downstream side intake pipe pressure Pm, 

[0056] Here, if making the average of the total amount of air flowing out from the intake pipe part 1 3' per unit time or 
the average of the amount of air taken in from the intake pipe part 13' to all combustion chambers 5 per unit time over 
the intake stroke of one cylinder the cylinder intake air flow mc (explained in detail below), since the cylinder filling air 

30 amount Mc is proportional to the downstream side intake pipe pressure Pm, the cylinder Intake air flow mc can also be * 
considered proportional to the downstream side intake pipe pressure Pm. From this, the equation (11 ) is obtained based 
on logic and experience. Note that the compliance parameter a in equation (1 1) is a proportional coefficient, while the 
compliance parameter b is a value relating to the amount of burnt gas remaining in the combustion chamber 5 when the 
exhaust valve is closed (explained below). 

35 [0057] Note that by having the compliance parameters a and b take two different values (for example, a1 , b1 and a2, 
b2) when the downstream side intake pipe pressure Pm is large and when it is small even if the engine speed etc. are 
the same, that is, by having the cylinder intake air flow mc shown by two equations like equation (11) (that is, the linear 
equation of downstream side intake pipe pressure Pm), it is learned that sometimes it is possible to find the cylinder 
intake air flow mc more accurately. This is believed to be related to the fact that, and in particular when both the intake 

^ valves 6 and exhaust valves 7 are open (that is, valve overlap) etc., the burnt gas flows back to the intake ports 7. That 
is, when there is valve overlap, when the downstream side intake pipe pressure Pm is the predetermined pressure or 
more, the higherthe downstream side intake pipe pressure Pm, the much less the backflowof the burnt gas, so compared 
with when the predetermined pressure or less, the value of a is made larger and the value of b is made smaller. 
[0058] Here, the cylinder intake air flow mc will be explained with reference to FIG. 8 for the case where the internal 

^5 combustion engine has four cylinders. Note that FIG. 8 shows on the abscissa the rotational angle of the crankshaft and 
the ordinate the amount of air actually flowing from the intake pipe part 13' to the combustion chamber 5 per unit time. 
As shown in FIG. 8, in a four-cylinder internal combustion engine, the intake valves 6 opens in the order of for example 
the #1 cylinder, #3 cylinder, #4 cylinder, and #2 cylinder and air flows from the intake pipe part 13' to the combustion 
chamber 5 of each cylinder in accordance with the amount of opening of the intake valves 6 corresponding to each 

50 cylinder. The change in the flow of the air flowing from the intake pipe part 13* to the combustion chamber 5 of each 
cylinder is as shown by the broken line in FIG. 8, while the flow of airf lowing from the intake pipe part 1 3' to the combustion 
chambers 5 of all cylinders Is as shown by the solid line in FIG. 8. Further, for example, the cylinder filling air amount 
Mc to the #1 cylinder corresponds to the part shown by the hatching in FIG. 8. 

[0059] As opposed to this, the average of the amount of air flowing from the intake pipe part 13' into the combustion 
55 chambers 5 of all cylinders shown by the solid line is the cylinder intake air flow mc and is shown by the one-dot chain 
line in the figure. Further, the cylinder intake air flow mc shown by the one-dot chain line multiplied with the time AT^so* 
taken for the crankshaft to rotate 180'' In the case of four cylinders (that is, in a four-stroke type internal combustion 
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engine, the angle 720** which the crankshaft rotates in one cycle divided by the number of cylinders) becomes the cylinder 
filling air amount Mc. Therefore, by multiplying the cylinder intake air flow mc calculated by the intake valve model M23 
with ATiBo-. 't is possible to calculate the cylinder filling air amount Mc {Mc=mc-AT^so»)* Further, by dividing this cylinder 
filling air amount Mc by the mass of the air occupying a volume conresponding to the amount of exhaust per cylinder in 
5 the state of 1 atm and 25'C, it is possible to calculate the cylinder air filling rate Kl. Note that, as clear from the above 
explanation, it is believed that if multiplying the value b in equation (1 1 ) with AT^qqc, the amount of bumt gas remaining 
in the combuslon chambers, at the time the exhaust valve 8 is closed, is obtained. 

[0060] Next, the case of using the intake air amount model M20 to actually calculate the cylinder filling air amount Mc 
will be explained. The cylinder filling air amount Mc is expressed by using the intake air amount model M20 to solve the 
10 equation (5), equation (7), equation (8), and equation (1 1). In this case, for processing at the ECU, it is necessary to 
make these equations discrete. If using the timet and calculation interval (discrete time) At to make equation (5), equation 
(7), equation (8), and equation (1 1 ) discrete, the following equation (12), equation (13), equation (14), and equation (15) 
are obtained. Note that, the downstream side intake pipe temperature Tm(t+At) is calculated by equation (16) from 
Pm/Tm(t+At) and Pm(t+At) calculated by equation (13)and equation (14), respectively: 



15 



20 



30 



35 



40 



mt(t)=M . At(0t(t)) .-^^L^^f^l (12) 

VR- Ta { Pa J 



— (t+At)=— (t)+Af — • (mt(t)-mc(t)) (13) 
Tm Tm Vm 



D 

Pm(t+At) = Pm(t)+At • s (mt(t) 'Ta-mcW • Tm(t)) (14) 

Vm 



mc(t) = a • Pm(t)-b (15) 



Tm(t+At)= y-*-^^^ (16) 



[0061] In the intake air amount model M20 mounted in this way, the throttle valve passage air flow mt(t) at the time t 
calculated by equation (12) of the throttle model M21 and the cylinder intake air flow mc(t) at the time t calculated by 
equation (15) of the intake valve model M23 are entered in equation (13) and equation (14) of the intake pipe model 
M22. Due to this, the downstream side intake pipe pressure Pm(t+At) and the downstream side intake pipe temperature 

45 Tm(t+At) at the time t+At are calculated. Next, the calculated Pm(t+At) are entered into equation (1 2) and equation (1 5) 
of the throttle model M21 and intake valve model M23. Due to this, the throttle valve passage air flow mt(t+At) and 
cylinder intake air flow mc(t+At) at the time t+At are calculated. Further, by repeating this calculation, the cylinder intake 
airflow mc at any timing t is calculated from the throttle valve opening degree Gt, atmospheric pressure Pa, and atmos- 
pheric temperature Ta, and the calculated cylinder intake airflow mc is multiplied with the time AT^goo so as to calculate 

50 the cylinder filling air amount Mc at any timing t. 

[0062] Note that, at the time of start of the internal combustion engine, that is, at the time t=0, the downstream side 
intake pipe pressure Pm is made equal to the atmospheric pressure (Pm(0)=Pa), the downstream side intake pipe 
temperature Tm is made equal to the atmospheric temperature (Tm(0)=Ta), and the calculations in the models M21 to 
M23 are started. 

55 [0063] Note that, in the intake air amount model M20, it is assumed that the atmospheric temperature Ta and atmos- 
pheric pressure Pa are constant, but it is also possible to make the values change along with time. For example, it is 
also possible to enterthe value detected at the time t by an atmospheric temperature sensor for detecting the atmospheric 
temperature as the atmospheric temperature Ta(t) and enter the value detected at the time t by an atmospheric pressure 
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10 



15 



20 



sensor for detecting the atmospheric pressure as the atnnospheric pressure Pa(t) Into the equation (12) and equation (14). 
[0064] However, when controlling an internal combustion engine, in particular when using a model to control an internal 
combustion engine, to calculate the parameters relating to control, sometimes the throttle vah/e downstream side intake 
pipe pressure Pmta and/or cylinder intake air flow mcta at the time of steady operation (or the cylinder air filling rate KIta 
at the time of steady operation able to be calculated from this) is necessary. Here, the value at the time of steady operation 
(Pmta, mcta, etc.) means the finally taken value when steadily operating the internal combustion engine in a certain 
state, that is, the value considered as the convergence value. These values are used in control of an Internal combustion 
engine mainly for avoiding complicated catoulatlon or reducing the amount of cateulation so as to reduce the control 
load or for improving the precision of the parameters calculated. Further, these values may be made the ones found 
using maps in the past. 

[0065] That is, for example, the practice Is to prepare In advance a map for finding that value using the throttle valve 
opening degree, engine speed, or other indicator of the operating state as arguments, store it in the ROM, and search 
through the map based on the operating state at that time to find the required value. However, to actually prepare such 
a map, a tremendous amount of time is required. That is, to prepare a map, it is necessary to actually measure the Pmta 
or mcta while successively changing the arguments. The work becomes tremendous. Further, there is the concern that 
an increase in the necessary maps or arguments will increase the map searching operation and increase the control load. 
[0066] Therefore, in the control system of an internal combustion engine of the present embodiment, when the Pmta 
and/or mcta (or KIta) is necessary, the method explained below is used to find it without using a map. Note that, as will 
be clear from the following explanation, this method utilizes the fact that at the time of steady operation, the throttle valve 
passage air flow mt and cylinder intake air flow mc match. 

[0067] That Is, the control system of an Internal combustion engine of the present embodiment provides as calculation 
equations of the throttle valve passage air flow mt the following equation (17) and equation (18) (that is, the equation 
(5) and equation (6). Below, "equation (17) etc.") 



25 



Pa 



(17) 
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[0068] Further, the control system of an internal combustion engine of the present embodiment provides as a calculation 
equation of the cylinder intake air flow mc the following equation (1 9) (that is, the equation (11)): 



45 



mc = a • Pm-b 



(19) 



50 



55 



[0069] Further, when the internal combustion engine is in steady operation, the throttle valve passage air flow mt and 
the cylinder Intake air flow mc match. Therefore, if finding the downstream side Intake pipe pressure Pm when the throttle 
valve passage airflow mt found from the equation (17) etc. and the cylinder intake airflow mc found from the equation 
(1 9) match, the downstream side Intake pipe pressure Pmta at the time of steady operation under the operating conditions 
at that time is found. Further, in the same way, by finding the cylinder intake airflow mc when the throttle valve passage 
air flow mt found from the equation (17) etc. and the cylinder intake airflow mc found from the equation (19) match, it 
is possible to find the cylinder intake air flow mcta at the time of steady operation under the operating conditions at that 
time (further, it is also possible to find the cylinder air filling rate KIta at the time of steady operation from this value). 
[0070] Further, finding the Pmta and mcta In this way is synonymous with finding the intersecting point EP between 
the curve mt expressed by the equation (17) etc. and the line mc expressed by the equation (19) as shown in FIG. 9. 
Here, when finding the intersecting point EP, if using equation (17) expressing the curve mt as it is to tiy to find the 
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intersecting point EP, the calculation becomes extremely complicated. Therefore, to simplify the calculation, it is also 
possible to approximate the equation (1 7) etc. by a plurality of linear equations of the downstream side intake pipe 
pressure Pm. That is, the curve mt is approximated by a plurality of lines. Specifically, for example, it is also possible to 
calculate the throttle valve passage airflow mt based on the equation (1 7) etc. at every constant interval of the downstream 
5 side intake pipe pressure Pm, find the points on the curve mt at constant intervals of the downstream side intake pipe 
pressure Pm, and find the lines connecting the adjoining two points as the approximated lines of the curve mt. Further, 
the linear equations expressing these approximated lines become the approximated linear equations, of the equation 
(17) etc. 

[0071] However, the approximation of the equation (17) etc. by linear equations is for easily finding the intersecting 
'0 point EP, so what is necessary here Is an approximated linear equation of the equation (17) etc. in the vicinity of the 
intersecting point EP. Therefore, it is also possible to find only this approximated linear equation. In this case, it is also 
possible to find in advance the cylinder intake air flow mc based on the equation (19) every constant interval of the 
downstream side intake pipe pressure Pm and find the point where the throttle valve passage air flow mt and the cylinder 
intake airflow mc invert in magnitude so as to specify the position of the intersecting point EP. 
IS [0072] More specifically, the approximated linear equation at the vicinity of the intersecting point EP (that is, the portion 
where the throttle valve passage air flow mt and cylinder intake air flow mc invert in magnitude), for example, is made 
the linear equation expressing the line nmt connecting the two points tj and tk on the curve mt expressed by the equation 
(1 7) etc. which are the points before and after the throttle valve passage air flow mt and the cylinder intake air flow mc 
invert in magnitude (see FIG. 10). 
20 [0073] Note that in the region where the downstream side intake pipe pressure Pm is the critical pressure Pc (that is, 
the pressure where even if the downstream side intake pipe pressure Pm becomes that pressure or less, the throttle 
valve passage flow mt does not increase further) or less, mt becomes a constant value, so even if not using the above- 
mentioned approximation, the intersecting point EP can be easily found. 

[0074] Further, when the compliance parameters a and b of equation (1 9) take two different values when the down- 
25 stream side intake pipe pressure Pm is large and when it is small (for example, a1 , b1 and a2, b2), that is, when as 
shown in FIG. 9 the cylinder intake airflow mc is shown by two lines connected by the connection point CP, when the 
connection point CP is in the vicinity of the intersecting point EP, by approximating the two lines by one line in the vicinity 
of the intersecting point EP, it is possible to simplify the calculations for finding the intersecting point EP and lighten the 
control load. 

30 [0075] Specifically, for example, as shown in FIG. 10, the two lines showing the cylinder intake air flow mc are ap- 
proximated by one line. That is, in this case, the cylinder intake airflow mc is shown by two equations expressed in the 
form of equation (1 9) (that is, two linear equations of downstream side intake pipe pressure Pm with different compliance 
parameters a and b), but these equations are approximated in the vicinity of the intersecting point EP by a linear equation 
expressing the line nmc connecting the points cj and ck, one of which is on one of the two lines mc expressed by the 

ss two equations and which are at the positions sandwiching the connection point CP and the intersecting point EP. 

[0076] In the example shown in FIG. 1 0, the curve mt showing the throttle valve passage air flow mt in the vicinity of 
the intersecting point EP is approximated by the line nmt, and the two lines expressing the cylinder intake air flow mc 
are approximated by a single line nmc. Due to this, the intersecting point nEP sought becomes slightly different from 
the intersecting point EP, but this intersecting point nEP can be simply found by calculation finding the intersecting point 

^ of the two lines nmt and nmc. That is, according to this method, it is possible to simply find the approximated values of 
the downstream side intake pipe pressure Pmta and cylinder intake air flow mcta at the time of steady operation. 
[0077] However, in the above-mentioned equation (1 7) etc., the throttle valve passage airflow mt is calculated using 
the intake pipe pressure at the upstream side of the throttle valve 1 8 (hereinafter referred to as the "upstream side intake 
pipe pressure") as the atmospheric pressure Pa. However, the actual upstream side intake pipe pressure usually becomes 

45 a pressure lower than atmospheric pressure during engine operation since there is pressure loss at the upstream side 
of the throttle valve in the engine intake system. In particular, in the conflguration shown in FIG. 1, the air cleaner 1 6 is 
provided at the upstream-most part in the engine intake system, so to more precisely calculate the throttle valve passage 
air flow mt, it is preferable to consider at least the pressure loss of the air cleaner 16. 

[0078] Therefore, in a control system of an internal combustion engine of another embodiment of the present invention, 
50 to more accurately calculate the throttle valve passage air flow mt, instead of the equafion (17) etc., it is also possible 
to provide the following equation (20) and equation (21 ) (hereinafter referred to as "equation (20) etc.") as the equations 
for calculation of the throttle valve passage airflow mt. In equation (20) etc., at the portion in the equation (1 7) etc. where 
the atmospheric pressure Pa is used, the upstream side intake pipe pressure Pac found considering at least the pressure 
loss of the air cleaner is used. 

55 
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[0079] By using the equation (20) etc. as the calculation equation of the throttle valve passage airflow mt, it is possible 
to more accurately find the downstream side intake pipe pressure Pmta and cylinder intake airflow mcta at the time of 
steady operation by the above-mentioned method. 

[0080] However, the upstream side Intake pipe pressure Pac may also be detected by providing a pressure sensor 
directly upstream of the throttle valve 18. but it Is also possible to calculate it without using a pressure sensor. That is, 
the difference between the atmospheric pressure Pa and the upstream side intake pipe pressure Pac can be expressed 
as in the following equation (22) by pBemouHi's theorum: 



25 



Pa-Pac = ipv='=k— 
2 p 



(22) 



30 



35 



[0081] Here, p is the atmospheric density, v is the velocity of the air passing through the air cleaner 1 6, Ga is the flow 
of the air passing through the air cleaner 1 6, and k is a proportional coefficient between v and Ga. If using the standard 
atmospheric density pO. the pressure connection coefficient ekpa for converting the standard atmospheric density pO to 
the current atmospheric density p, and the temperature correction coefficient ektha, equation (22) can be replaced by 
equation (23). Further, equation (23) can be replaced by equation (24) using the function f(Ga) having just the flow Ga 
as a variable. 



Pa-Pac = — • Ga^ 

p 0 ekpa • 



ektha 



40 
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[0082] Equation (24) can be modified to equation (25) showing the upstream side intake pipe pressure Pac. In equation 
(25), the flow Ga can be detected by an air flow meter provided at the immediately downstream side of the air cleaner 
1 6 when such an air flow meter Is provided. Further, the pressure correction coefficient ekpa can be set by the detected 
atmospheric pressure Pa, while the temperature coaectlon coefficient ektha can be set by the detected atmospheric 
temperature Ta. 



Pac = Pa- 



f(Ga) 



ekpa * ektha 



(25) 



[0083] Further, in equation (25), the flow Ga of the air passing through the air cleaner 1 6 can be considered to be the 
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throttle valve passage air flow mt Equation (25) can be modified to equation (26): 

Pac = Pa (26) 

ekpa • ektha 

[0084] However, since the current upstream side intake pipe pressure Pac is necessary for calculating the current 
10 throttle valve passage air flow mt based on equation (20) etc., to calculate the current upstream side intake pipe pressure 
Pac based on equation (26), it is necessary to use as the throttle valve passage air flow mt, the previous throttle valve 
passage air flow mt, that is, throttle valve passage air flow mt of one discrete time before. In this respect, repeated 
calculation can improve the precision of the calculated upstream side intake pipe pressure Pac, but to avoid an increase 
in control load, it is also possible to use the upstream side intake pipe pressure Pac found based on the previously found 
15 throttle valve passage air flow mt as the cun^ent (present) upstream side intake pipe pressure Pac. 

[0085] Further, it is also possible to use the following method to find the downstream side intake pipe pressure Pmta 
and the cylinder intake air flow mcta at the time of steady operation when considering at least the pressure loss of the 
air cleaner 16. That is, with this method, in the above method of approximating the equation (17) etc. at least in the 
vicinity of the intersecting point EP by a linear equation, finding the intersecting point between the approximated line 
20 shown by that approximated linear equation and the line expressed by the equation (1 9) (or its approximated line), and 
finding the downstream side intake pipe pressure Pmta and cylinder intake airflow mcta at the time of steady operation, 
the approximated linear equation of the equation (17) etc. (or approximated line expressed by the approximated linear 
equation) is corrected using the upstream side intake pipe pressure Pac. 

[0086] That Is, in the above method, the approximated line of the curve mt expressed by the equation (17) etc., as 
25 shown in FIG. 10, was found as the line nmt connecting the two points tj and tk on the curve mt before and after the 
throttle valve passage air flow mt and cylinder intake air flow mc invert in magnitude, but with this method, the values 
of the downstream side intake pipe pressure and throttle valve passage air flow showing the coordinates of the two 

points tj and tk are multiplied with Pac/Pa and the line connecting the two points shown by the new coordinates (approx- 
imated line after correction) is found (the linear equation expressing this line becomes a corrected approximated linear 
30 equation). 

[0087] Further, by finding the intersecting point between this approximated line after connection and the line expressed 
by the equation (19) (or its approximated line), the downstream side intake pipe pressure Pmta and cylinder intake air 
flow mcta at the time of steady operation when considering at least the pressure loss of the air cleaner 1 6 are found. 
[0088] Next, another embodiment of the present invention will be explained with reference to FIG. 11. FIG. 1 1 is a 
35 schematic view of an example of the case of applying a control system of an internal combustion engine of the present 
invention to a cylinder injection type spark ignition intemal combustion engine different from FIG. 

1. The configuration shown in FIG. 1 1 is basically the same as the configuration shown in FIG. 1. Explanations of 
common parts are in principle omitted. 

40 

[0089] Compared with the configuration shown in FIG. 1 , the configuration shown in FIG. 1 1 differs in that the exhaust 
passage (exhaust port, exhaust pipe, etc.) and intake passage (intake port and intake pipe) are connected to each other 
through an exhaust gas recirculation passage (hereinafter refen^d to as "EGR passage") 21 and a control valve 22 for 
adjusting the flow of the exhaust gas passing through the exhaust gas recirculation passage 21 (hereinafter referred to 
45 as an "EGR control valve") is arranged in this exhaust gas recirculation passage 21 . That is, in the present embodiment, 
sometimes exhaust gas recirculation for making part of the exhaust gas discharged to the exhaust passage flow into 
the intake passage (hereinafter referred to as "EGR") is perfomied. 

[0090] Further, the configuration shown in FIG. 1 1 differs from the configuration shown in FIG. 1 In the point that it is 
provided with a variable valve timing mechanism 23 for changing the operating timing of the intake valves 6. Note that 

50 the EGR control valve 22 and variable valve timing mechanism 23 are both controlled by the ECU 31 . 

[0091] Further, in the present embodiment as well, a model is constmcted for the configuration shown in FIG. 1 1. In 
the same way as the case of the above-mentioned other embodiment, the model is used for control of the intemal 
combustion engine. Further, in the present embodiment as well, in the same way as the above-mentioned other em- 
bodiment, when the downstream side intake pipe pressure Pmta and/or cylinder intake air flow mcta at the time of steady 

55 operation (or the cylinder air filling rate KIta at the time of steady operation able to be calculated from this) is necessary, 
the fact that the throttle valve passage air flow mt and the cylinder intake air flow mc match at the time of steady operation 
is utilized and these values are found by calculation. 

[0092] However, in the present embodiment, EGR is sometimes perfomied. Further, sometimes the operating timing 
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of the intake valves 6 (hereinafter referred to as simply as "valve timing") is changed. For this reason, the control system 
of an internal combustion engine of the present embodiment is provided with the following equation (27) instead of the 
equation (1 9) as the calculation equation of the cylinder intake air flow mc used for calculation of the Pmta and/or mcta. 
[0093] That is. in the present embodiment, the downstream side intake pipe pressure Pm when the throttle valve 
5 passage air flow mt found from the equation (1 7) etc. and cylinder intake air flow mc found from the following equation 
(27) match is found as the Pmta, and the cylinder intake air flow mc at that time is found as the mcta. Alternatively, when 
considering at least the pressure loss due to the air cleaner 16, the downstream side intake pipe pressure Pm when the 
throttle valve passage air flow mt found from the equation (20) etc. and the cylinder intake air flow mc found from the 
following equation (27) match is found as the Pmta, and the cylinder intake air flow mc at that time is found as the mcta. 

10 

mc = e • Pm+g (27) 

[0094] Equation (27) is an equation obtained since even if EGR is perfomned and/or the valve timing is changed, the 
cylinder intake air flow mc changes substantially linearly based on the downstream side intake pipe pressure Pm. Here, 
e and a are compliance parameters different from the compliance parameters a and b in equation (19) (or equation 
(11)), that is, are compliance parameters determined based on at least the engine speed NE, EGR control valve opening 
degree STP, and valve timing VT. Further, by making the compliance parameters e and g take different values for each 
predetennined range of the downstream side intake pipe pressure Pm even if the engine speed NE, EGR control valve 
^ opening degree STP, valve timing VT, or other operating conditions are the same, that is, by expressing the cylinder 
intake airflow mc by a plurality of equations like the equation (27) (that is, linear equation of downstream side intake 
pipe pressure Pm), it is learned that sometimes the cylinder intake air flow mc can be found more accurately. 
[0095] The compliance parameters e and may be found by preparing in advance a map using the engine speed N E, 
EGR control valve opening degree STP, and valve timing VT as arguments and searching through the map based on 
the operating conditions at that time in accordance with need, but it is also possible to use the method explained below 
to estimate the necessary compliance parameters e and 3 and therefore sharply cut the manhours for making the map. 
Further, if using this method in accordance with need to estimate the compliance parameters e and 3, it is possible to 
reduce the number of maps stored and lighten the control load for map searching. 

[0096] That is, with this method, for the case of each engine speed NE, only the compliance parameters exn and gxn 
30 when the EGR control valve opening degree STP is a certain single EGR control valve opening degree STPx and when 
making the valve timing VT each valve timing VTn and the compliance parameters emx and gmx when the valve timing 
VT is a certain single valve timing VTx and when making the EGR control valve opening degree STP each EGR control 
valve opening degree STPm are found and these are used to estimate compliance parameters emn and gmn at the 
time of any other EGR control valve opening degree STPm and any other valve timing VTn. This method utilizes the 
^ fact that when the engine speed NE is constant, the amount of the EGR gas sucked into the cylinder is substantially 
detennined by the EGR control valve opening degree STP and downstream side intake pipe pressure Pm. 
[0097] Below, this will be explained more specifically. For example, when the engine speed NEisNEI , if the compliance 
parameters when the EGR control valve opening degree STP is the closed state STPO and the valve timing VT is the 
reference timing VTO (that is, advance =0) are designated as eOO. gOO, the cylinder intake airflow mcOO at that time can 
^ be expressed by the following equation (26). 



mcOO = eOO • Pm+gOO (28) 

[0098] Similarly, if the compliance parameters when the EGR control valve opening degree STP is STP1 and the 
valve timing VT is the reference timing VTO (that is, advance=0) are designated as elO, glO, the cylinder intake air flow 
mclO at that time can be expressed by the following equation (29). 



mclO = elO • Pm+glO (29) 

[0099] Further, from these equation (28) and equation (29), the flow mcegri 0 of the EGR gas flowing into the cylinder 
(hereinafter refen-ed to as "cylinder intake EGR flow") when the EGR control valve opening degree STP is STP1 and 
the valve timing VT is the reference timing VTO (that is, advance =0) can be expressed by the following equation (30). 
Here E and G are coefficients expressing the calculated values of the corresponding compliance parameters. 
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racegrlO =nicOO-mclO 

= (eOO-elO) • Pm +(gOO-glO) = E • Pm+G 



(30) 



[0100] If illustrating these equation (28) to equation (30), for example, the result becomes FIG. 12. In the example of 
FIG. 12, the compliance parameters el 0 and g1 0 are assumed to take different values when the downstream side intake 
pipe pressure Pm is the predetermined pressure Pm1 or more and when it is less than the predetemriined pressure Pml. 
As a result, the coefficients E and G are assumed to take different values when the downstream side intake pipe pressure 
Pm is the predetemriined pressure Pm1 or more and when it is less than the predetermined pressure Pnil. Further, in 
this example, when the downstream side intake pipe pressure Pm is less than the predetermined pressure Pm1, the 
compliance parameters eOO and el 0 are assumed to be substantially equal. 

[01 01 ] Further, in the same way as eq uatlon (28) and equation (29), when the engine speed N E is N El , if the compliance 
parameters when the EGR control valve opening degree STP is the closed state STPO and the valve timing VT is VT1 
are designated as eOl , gOI , the cylinder intake airflow mcOl at that time can be expressed by the following equation (31 ). 



[0102] Further, here, when the engine speed NE is constant, the amount of the EGR gas taken into the cylinder is 
substantially detemiined by the EGR control valve opening degree STP and the downstream side intake pipe pressure 
Pm. If considering this, the cylinder intake EGR flow mcegrl 1 when the EGR control valve opening degree STP is STP1 
and the valve timing VT is VT1 is substantially equal to the above mcegrl 0 and can be expressed by the above equation 



[0103] Further, from this, the cylinder intake airflow mcl 1 when the EGR control valve opening degree STP is STP1 
and the valve timing VT is VT1 can be expressed as in the following equation (32) from equation (30) and equation (31 ). 



[0104] That is, the compliance parameters e1 1 and g1 1 when the EGR control valve opening degree STP is STP1 
and the valve timing VT is VT1 are expressed as in the following equation (33). That is, the compliance parameters e1 1 
and g1 1 when the EGR control valve opening degree STP is STP1 and the valve timing VT is VT1 can be estimated 
from the compliance parameters eOO and gOO when the EGR control valve opening degree STP is STPO and the valve 
timing VT is VTO, the compliance parameters e10 and glO when the EGR control valve opening degree STP is STP1 
and the valve timing VT is VTO, and the compliance parameters e01 and gOI when the EGR control valve opening 
degree STP is STPO and the valve timing VT is VT1 . 



[0105] These equation (30), equation (31), and equation (32), if illustrated, becomes as shown for example In FIG. 
13. In the example of FIG. 13, the compliance parameters eOl and gOI are assumed to take different values when the 
downstream side intake pipe pressure Pm is the predetennined pressure Pm2 or more and when it Is less than the 
predetemriined pressure Pm2. 

[01 06] Note that, above, to simplify the explanation, the explanation was given of the case of estimating the unknown 
compliance parameters e1 1 and g1 1 based on the case where the EGR control valve opening degree STP was in the 
closed state STPO, but the present invention is not limited to this. IHowever, when the EGR control valve opening degree 
STP is in the closed state STPO, compared with other cases, it is possible to find the compliance parameters e and ^ 
more precisely, so by using the case where the EGR control valve opening degree STP is the closed state STPO as the 
standard, it becomes possible as a result to precisely estimate the unknown compliance parameters el 1 and g1 1 . 
[0107] Further, as clear from the above explanation, according to this method, if finding the compliance parameters 
exn and gxn when the EGR control valve opening degree STP is a certain EGR control valve opening degree STPx and 



nic01 = e01 • Pm+gOl 



(31) 



(30). 



mcll = (e01-eOO + elO) • Pm +(g01-gOO + glO) 
= (e01-E) •Pm + (g01-G) 



(32) 




(33) 
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the valve timing VT is each valve timing VTn and compliance parameters emx and gmx when the valve timing VT is a 
certain valve timing VTx and the EGR control valve opening degree STP is each EGR control valve opening degree 
STPm for the case of each engine speed NE, it is possible to use these to estimate the compliance parameters emn 
and gmn at the time of any other EGR control valve opening degree STPm and any valve timing VTn. Further, from this, 
it is possible to greatly reduce the manhours for map making. 

[0108] However, as in the case shown in FIG. 13, when both the cylinder intake EGR flow mcegrlO and the cylinder 
intake air flow mcOl are shown by two lines connected by connection points, the cylinder intake airflow mc1 1 estimated 
based on these becomes shown by three lines connected by two connection points. If the cylinder intake air flow is 
shown by three lines in this way, compared with when shown by two lines, the processing for calculation when finding 
the intersecting point with the curve etc. showing the throttle valve passage air flow mt for finding the Pmta and/or mcta 
becomes extremely troublesome. 

[0109] Therefore, to lighten the control load, it is also possible to approximate the three lines showing the cylinder 
intake air flow by two lines by the method explained below. That is, this method approximates the three lines expressing 
the estimated cylinder intake airflow mc1 1 by two lines using as a reference point, of the two connection points connecting 
them, the connection point RP with the same Pm coordinate as the connection point of the two lines expressing the 
cylinder intake air flow mcOl forming the basis for the estimation. That is, an equation expressing the two approximated 
lines connected by the connection point RP is found. What is expressed by these two lines is the approximated cylinder 
intake airflow mc'l 1 approximating the cylinder intake airflow mc1 1 . Below, this will be explained more specifically with 
reference to FIG. 14 and FIG. 15. 

[0110] As shown in FIG. 14 and FIG. 15, when the cylinder intake EGRflowmcegrlOis shown by two lines connected 

by a connection point, In the equation (30), the coefficients E and G take different values when the downstream side 
intake pipe pressure Pm is the predetemriined pressure Pml or more and when it is less than the predetermined pressure 
Pml. In this case, if the cylinder intake EGR flow mcegrlO when the downstream side intake pipe pressure Pm is less 
than the predetemriined pressure Pml is mcegrllG, the coefficients E and G when the downstream side intake pipe 
pressure Pm is less than the predetermined pressure Pml are El and Gl, the cylinder intake EGR flow mcegrlO when 
the downstream side intake pipe pressure Pm is the predetermined pressure Pm1 or more is mcegrhIO, and the coef- 
ficients E and G when the downstream side intake pipe pressure Pm is the predetennined pressure Pm1 or more are 
Eh and Gh, the equation (30) can be expressed as in the following equation (34). 



rocegrllO =E1 • Pm+Gl, Pm<Pml 1 

racegrhlO = Eh • Pm + Gh, Pm ^ Pml J ^ ^ ^ 

[0111] Similarly, as shown in FIG. 14 and FIG. 15, when the cylinder intake air flow mcOl is shown by two lines 
connected by a connection point, the compliance parameters eOl and g01 in the equation (31) take different values 
when the downstream side intake pipe pressure Pm is a predetennlned pressure Pm2 or more and when it is less than 
a predetermined pressure Pm2. In this case, if the cylinder intake air flow mcOl when the downstream side intake pipe 
pressure Pm is less than the predetennined pressure Pm2 is mcllO, the compliance parameters eOl and gOI when the 
downstream side intake pipe pressure Pm is less than the predetemiined pressure Pm2 are el 01 and g101 , the cylinder 
intake air flow mcOl when the downstream side intake pipe pressure Pm is the predetennined pressure Pm2 or more 
is mchOI, and the compliance parameters eOl and g01 when the downstream side intake pipe pressure Pm is the 
predetermined pressure Pm2 or more are eh01 and gh01 , the equation (31 ) can be expressed by the following equation 
(35): 



racl01 = el01 • Pm+glOl, Pm<Pm2 \ 
mchOl = ehOl • Pm +gh01, Pm ^ Pm2 J 

[0112] Further, with this method, the cylinder intake air flow mcl 1 is approximated by a line of the slant (el 01 -El) 
when the downstream side intake pipe pressure Pm is less than a predetennined pressure Pm2 and is approximated 
by a line of the slant (ehOI -Eh) when the downstream side intake pipe pressure Pm is the predetennined pressure Pm2 
or more. Further, these two approximated lines are connected by the connection point RP. 

[0113] The equation expressing this approximated line, that is, the equation expressing the approximated cylinder 
intake airflow mc'1 1 approximating the cylinder intake air f low mc1 1, becomes a different equation when the downstream 
side intake pipe pressure Pm is the predetermined pressure Pm2 or more and when It is less than the predetermined 
pressure Pm2 and is used selectively depending on the relative magnitude of the predetemiined pressures Pml and Pm2. 
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[011 4] The equation expressing the approximated cylinder intake airflow mc'1 1 found by this method, when Pml >Pm2 
as shown in FIG. 14, can be expressed by the following equation (36) when the approximated cylinder intake air flow 
mc'1 1 when the downstream side intake pipe pressure Pm is less than the predetemiined pressure Pm2 Is mc'1 1 1 and 
the approximated cylinder intake airflow mcM 1 when the downstream side intake pipe pressure Pm is the predetermined 
pressure Pm2 or more is mc'hl 1 : 



mc Ill = (el01-El) • Pra + (gl01-Gl) 

= eplall • Pm +gplall, Pm <Pm2 
mc'hll = (eh01-Eh) •Pm + (Eh-El) • Pm2+(gh01-Gl) 



= (eh01-Eh) 
= ephall • Pm +gphall, 



p„^(Eh-El)(gh01-glOl) 

elOl-ehOl 



> (36) 



Pm ^Pm2 



[0115] Here, epiall, gplal 1 . epha1 1 , and gphal 1 are coefficients obtained by rewriting the con^esponding parts in 
the equations and are approximated compliance parameters. Further, in this case, the coordinates of the connection 
point RP in FIG. 14can be expressed using the predetermined pressure Pm2 as (Pm2, (eh01-EI)'Pm2+(gh01-GI)). 
[01 1 6] On the other hand, as shown in FIG. 1 5, when Pm1 <Pm2, the equation expressing the approximated cylinder 
intake air flow mc'1 1 can be expressed by the following equation (37). 



mc Ill=(el01-El) • Pm+(eh01-el01+El-Eh)Pm2+(gh01-Gh) 

= (sl01-El) • P..'"'-"°'^^'-°-""""-'l°».(gh01-Gh) 

elOl-ehOl ^^3^^ 

= eplbll • Pm+gplbll, PmC:Pm2 
mchll = (eh01-Eh) • Pm+(gh01-Gh) 

= ephbl 1 • PnH-gphbl 1, PmS Pm2 



Here, epibi 1 , gplbl 1 , ephbl 1 , and gphbl 1 are coefficients obtained by rewriting the corresponding parts of the equations 
and are approximated compliance parameters. Further, in this case, the coordinates of the connection point RP in FIG. 
15 can be expressed using the predetermined pressure Pm2 as (Pm2, (e101-Eh) Pm2+(g101-Gh)). 
[0117] Further, as clear from FIG. 14 and FIG. 15, when using this method to find the approximated cylinder intake 
airflow mc'1 1 , when Pm1>Pm2, if the downstream side intake pipe pressure Pm is less than the predetermined pressure 
Pm2, the approximated cylinder intake air flow mc'1 1 matches with the cylinder intake air flow mcl 1 , while when 
Pml <Pm2, if downstream side intake pipe pressure Pm Is the predetermined pressure Pm2 or more, the approximated 
cylinder intake air flow mc'1 1 matches with the cylinder intake air flow mcl 1 . Note that, when Pml =Pm2, as the cylinder 
Intake air flow mc1 1 is originally shown by two lines, it is not necessary to use the above-mentioned method to find the 
approximated cylinder intake air flow me'1 1 . 

[0118] Further, it Is learned that even if using the method explained above to find the approximated cylinder intake air 
flow mc'1 1 and using this as the basis to find the Pmta or mcta, the effect on the calculation precision is relatively small. 
This is because when the engine speed NE is from the low speed to medium speed, the tendency is that Pm1=Pm2, 
while when the engine speed NE is the high speed, the tendency is that ElaEh. 

[0119] Note that, only naturally, even in the case of EGR explained later, it is possible to linearly approximate the 
curve showing the throttle valve passage airflow mt by any one of the methods explained above to find the Pmta or mcta. 
[0120] Further, in the configuration shown in FIG. 1 1, the variable valve timing mechanism 23 was provided only at 
the intake valve 6 side, but the present invention is not limited to this. That is, for example, a variable valve timing 
mechanism may also be provided at only the exhaust valve 8 side or may be provided at both of the intake valve 6 side 
and the exhaust valve 8 side. 

[0121] Further, the configuration shown in FIG. 1 1 has a variable valve timing mechanism 23 as an example of a 
variable intake apparatus, but the present invention can also be applied to cases of other variable intake apparatuses, 
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for example swirl control valves. That is, for example, for estimation of the compliance parameters e and g of equation 
(27). in the same way as the above-mentioned method, it is possible to estimate compliance parameters emn and gmn 
at the time of any EG R control valve opening degree STPm and any swirl control valve state SCn, for each engine speed 
NE, from the compliance parameters eyn and gyn when the EGR control valve opening degree STP is a certain EGR 
control valve opening degree STPy and the swirl control valve is each state SCn and the compliance parameters emy 
and gmy when the swirl control valve Is a certain state SCy and the EGR control vah/e opening degree STP is each 
EGR control valve opening degree STPm. 

[01 22] Note that the present invention was explained in detail based on specific embodiments, but a person skilled In 
the art could make various modifications, con-ections, etc. without departing from the claims and concept of the present 
invention. 



1. A control system of an internal combustion engine provided with: 

a throttle valve passage air flow calculation equation by which a throttle vafve passage air flow is expressed as 
a function of a downstream side Intake pipe pressure at the downstream side of a throttle valve and 
a cylinder intake air flow calculation equation by which a cylinder intake air flow is expressed as a function of 
said downstream side intake pipe pressure, 

said downstream side intake pipe pressure when the throttle valve passage air flow found from said throttle 
valve passage air flow calculation equation and the cylinder intake air flow found from said cylinder intake air 
flow calculation equation match being calculated as the downstream side intake pipe pressure at the time of 
steady operation under the operating conditions at that time. 

2. A control system of an internal combustion engine provided with: 

a throttle valve passage airflow calculation equation by which a throttle valve passage air flow is expressed as 
a function of a downstream side intake pipe pressure at the downstream side of a throttle valve and 
a cylinder intake air flow calculation equation by which a cylinder intake air flow is expressed as a function of 
said downstream side intake pipe pressure, 

said cylinder intake air flow when the throttle valve passage air flow found from said throttle valve passage air 
flow calculation equation and the cylinder intake air flow found from said cylinder intake air flow calculation 
equation match being calculated as the cylinder intake air flow at the time of steady operation under the operating 
conditions at that time. 

3. A control system of an internal combustion engine as set forth In claim 1 , wherein said cylinder intake air flow when 
a throttle valve passage air flow found from said throttle valve passage air flow calculation equation and a cylinder 
intake air flow found from said cylinder intake air flow calculation equation match is calculated as the cylinder intake 
air flow at the time of steady operation under the operating conditions at that time. 

4. A control system of an internal combustion engine as set forth in claim 1 , wherein 

said throttle valve passage airflow calculation equation is expressed as the following equation (1) where mt is a 
throttle valve passage air flow, p. is a flow coefficient at the throttle valve, At is a cross-sectional area of the opening 
of the throttle valve, Pa is an atmospheric pressure, Ta is an atmospheric temperature, R is a gas constant, Pm is 
said downstream side intake pipe pressure, and 0(Pm/Pa) is a coefficient determined in accordance with the value 
of Pm/Pa, and 

said cylinder intake air flow calculation equation is expressed as the following equation (2) where mc is a cylinder 
intake air flow and a and b are compliance parameters detemiined based on at least the engine speed: 



Claims 




(1) 
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mc = a • Pm-b (2) 



A control system of an internal combustion engine as set forth in claim 2, wherein 

said throttle valve passage air flow calculation equation is expressed as the following equation (1) where mt is a 
throttle valve passage airflow, is a flow coefficient at the throttle valve, At is a cross-sectional area of the opening 
of the throttle valve, Pa is an atmospheric pressure, Ta Is an atmospheric temperature, R is a gas constant, Pm is 
said downstream side intake pipe pressure, and <I>(Pm/Pa) Is a coefficient detenmined In accordance with the value 
of Pm/Pa, and 

said cylinder intake air flow calculation equation is expressed as the following equation (2) where mc Is a cylinder 
intake air flow and a and b are compliance parameters determined based on at least the engine speed: 



mt = ^ .At- Pi=.$f^l (1) 
VR-Ta \PaJ 



mc = a • Pm-b (2) 



A control system of an internal combustion engine as set forth in claim 1 , wherein said internal combustion engine 
has an exhaust gas recirculation passage for making at least part of the exhaust gas discharged into the exhaust 
passage flow into the Intake passage and an EGR control valve for adjusting the flow of the exhaust gas passing 
through said exhaust gas recirculation passage, 

said throttle valve passage air flow calculation equation is expressed as the following equation (3) wherein mt is a 
throttle valve passage air flow, |ji is a flow coefficient at the throttle valve, At is a cross-sectional area of the opening 
of the throttle valve. Pa Is an atmospheric pressure, Ta is an atmospheric temperature, R Is a gas constant, Pm is 
said downstream side intake pipe pressure, and <t>(Pm/Pa) is a coefficient detemiined in accordance with the value 
of Pm/Pa, and 

said cylinder intake air flow calculation equation Is expressed as the following equation (4) where mc is a cylinder 
intake air flow, and e and g are compliance parameters detemiined based on at least the engine speed and the 
opening degree of said EGR control valve: 



„,t = M -Af-^-O^I^) (3) 



Pa 
VR • Ta 



inc = e • Pm + g (4) 



A control system of an internal combustion engine as set forth in claim 2, wherein said internal combustion engine 
has an exhaust gas recirculation passage for making at least part of the exhaust gas discharged into the exhaust 
passage flow into the intake passage and an EGR control valve for adjusting the flow of the exhaust gas passing 
through said exhaust gas recirculation passage, 

said throttle valve passage air flow calculation equation is expressed as the following equation (3) where mt is a 
throttle valve passage air flow, \x is a flow coefficient at the throttle valve, At is a cross-sectional area of the opening 
of the throttle valve, Pa is an atmospheric pressure, Ta Is an atmospheric temperature, R is a gas constant, Pm is 
said downstream side intake pipe pressure, and 0(Pm/Pa) is a coefficient detemiined in accordance with the value 
of Pm/Pa, and 

said cylinder intake air flow calculation equation is expressed as the following equation (4) where mc is a cylinder 
intake air flow, and e and g are compliance parameters detennined based on at least the engine speed and the 
opening degree of said EGR control valve: 
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mt=fi • At •• 



VR • Ta 



Pa 




(3) 



5 



mc = e • Pm + g 



(4) 



8. A control system of an interna! combustion engine as set forth in claim 6, wherein said internal combustion engine 
further has a variable valve timing mechanism for changing an operating timing of a valve provided in each cylinder 
and, based on said compliance parameters e and g when said operating timing is a first valve timing and said EGR 
control valve is at a first opening degree, said compliance parameters e and g when said operating timing is said 
first valve timing and said EGR control valve is at a second opening degree, and said compliance parameters e and 

^5 2 when said operating timing is a second valve timing and said EGR control valve is at a first opening degree, said 

compliance parameters e and q when said operating timing is said second valve timing and said EGR control valve 
is at said second opening degree are estimated. 

9. A control system of an internal combustion engine as set forth in claim 7, wherein said internal combustion engine 
^ further has a variable valve timing mechanism for changing an operating timing of a valve provided in each cylinder 

and, based on said compliance parameters e and g when said operating timing is a first valve timing and said EGR 
control valve Is at a first opening degree, said compliance parameters e and g when said operating timing is said 
first valve timing and said EGR control valve is at a second opening degree, and said compliance parameters e and 
2 when said operating timing is a second valve timing and said EGR control valve is at a first opening degree, said 
compliance parameters e and ^ when said operating timing is said second valve timing and said EGR control valve 
is at said second opening degree are estimated. 

10. A control system of an internal combustion engine as set forth in claim 8, wherein when said compliance parameters 
e and when said operating timing is said second valve timing and said EGR control valve is at a first opening 

^ degree, respectively, take two different values when said throttle valve downstream side intake pipe pressure is 

larger than and smaller than a first pressure and 

said compliance parameters e and g when said operating timing is said second valve timing and said EGR control 
valve is at said second opening degree are estimated to take three or more different values in accordance with said 
throttle valve downstream side intake pipe pressure, 

^ based on said compliance parameters e and g when said operating timing is a first valve timing and said EGR control 

valve is at a first opening degree, said compliance parameters e and 2 when said operating timing is said first valve 
timing and said EGR control valve is at a second opening degree, and said compliance parameters e and g when 
said operating timing is a second valve timing and said EGR control valve is at a first opening degree, approximated 
compliance parameters ep and gp designed to take two values differing when said throttle valve downstream side 

^ intake pipe pressure is larger and smaller than a first pressure are calculated and these are made said compliance 

parameters e and £ when said operating timing is said second valve timing and said EGR control valve is at said 
second opening degree. 

11 . A control system of an internal combustion engine as set forth in claim 9, wherein when said compliance parameters 
^ e and c| when said operating timing is said second valve timing and said EGR control valve is at a first opening 

degree respectively take two values different when said throttle valve downstream side intake pipe pressure is larger 
than and smaller than a first pressure and 

said compliance parameters e and g when said operating timing is said second valve timing and said EGR control 
valve is at said second opening degree are estimated to take three or more different values in accordance with said 

^ throttle valve downstream side intake pipe pressure, 

based on said compliance parameters e and £ when said operating timing is a first valve timing and said EG R control 
valve is at a first opening degree, said compliance parameters e and g when said operating timing is said first valve 
timing and said EGR control valve is at a second opening degree, and said compliance parameters e and £ when 
said operating timing is a second valve timing and said EGR control valve is at a first opening degree, approximated 

^ compliance parameters ep and gp designed to take two different values when said throttle valve downstream side 

intake pipe pressure is largerthan and smallerthan a first pressure are calculated and these are made said compliance 
parameters e and £ when said operating timing is said second valve timing and said EGR control valve is at said 
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second opening degree. 

1 2. A control system of an internal combustion engine as set forth in claim 8, wherein the case where said EGR control 
valve is at said first opening degree is the case where said EGR control valve is closed. 

13. A control system of an internal combustion engine as set forth in claim 9, wherein the case where said EGR control 
valve is at said first opening degree is the case where said EGR control valve is closed. 

14. A control system of an internal combustion engine as set forth in claim 4, wherein at a portion where the throttle 
valve passage air flow mt and cylinder intake air flow mc invert in magnitude, said throttle valve passage air flow 
calculation equation used is an approximation equation expressed by a linear equation of the downstream side 
intake pipe pressure Pm. 

15. A control system of an internal combustion engine as set forth in claim 5, wherein at a portion where the throttle 
valve passage air flow mt and cylinder intake air flow mc invert in magnitude, said throttle valve passage air flow 
calculation equation used is an approximation equation expressed by a linear equation of the downstream side 
intake pipe pressure Pm. 

16. A control system of an internal combustion engine as set forth in claim 14, wherein said approximation equation is 
made a linear equation expressing a line connecting two points on a curve expressed by said throttle valve passage 
air flow calculation equation and before and after the point where the throttle valve passage air flow mt and cylinder 
intake air flow mc invert in magnitude. 

17. A control system of an internal combustion engine as set forth in claim 15, wherein said approximation equation is 
made a linear equation expressing a line connecting two points on a curve expressed by said throttle valve passage 
air flow calculation equation and before and after the point where the throttle valve passage airflow mt and cylinder 
intake air flow mc invert in magnitude. 

18. A control system of an intemal combustion engine as set forth in claim 4, wherein instead of said atmospheric 
pressure Pa, a throttle valve upstream side intake pipe pressure Pac found considering at least a pressure loss of 
an air cleaner is used. 

19. A control system of an intemal combustion engine as set forth in claim 5, wherein instead of said atmospheric 
pressure Pa, a throttle valve upstream side intake pipe pressure Pac found considering at least a pressure loss of 
an air cleaner is used. 

20. A control system of an internal combustion engine as set forth in claim 14, wherein 

a throttle valve upstream side Intake pipe pressure Pac found considering at least a pressure loss of an air cleaner 
is found based on the previously found throttle valve passage air flow, and 

said approximation equation is made a linear equation expressing a line connecting two points shown by coordinates 
obtained by multiplying Pac/Pa with values of the downstream side intake pipe pressure and the throttle valve 
passage air flow showing coordinates of two points on a curve expressed by said throttle valve passage air flow 
calculation equation and before and after a point where a throttle valve passage airflow mt and cylinder intake air 
flow mc invert in magnitude. 
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Fig. 5 
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Fig. 9 
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Fig. 11 
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I . . • engine body 

5. . .combustion chamber 
6. . .intake valve 
7. . .intake port 
8. . .exhaust valve 
9... exhaust port 

II . . . fuel injector 
13. . . intake pipe 

18 ... throttle valve 

22...EGR control valve 

23. . .variable valve timing mechanism 
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